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Abstract 

We analyze the signal at the Tevatron for the doubly charged Higgs super- 
field of supersymmetric left-right models with gauge mediated supersymmetry 
breaking. Due to the presence of a new coupling, the lighter stau is usually 
less massive than the lightest neutralino and is frequently the NLSP. The 
fermionic part of the doubly charged Higgs field decays dominantly into two 
r leptons plus missing energy (gravitino). We find that the inclusive two t- 
jets and three r-jets plus missing energy signatures could be observable at 
Run II of the Tevatron. We also find that the distribution in angle between 
the two highest Er r-jets when they come from same sign r leptons can be 
used to distinguish this model from the others. 
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I. INTRODUCTION 



Due to its many attractive features, supersymmetry (SUSY) has become a primary focus 
of experimental investigations. Most of the experimental searches for SUSY so far have 
been performed in the context of the Minimal Supersymmetric Standard Model (MSSM) 
IJ. It is important, however, to determine what the signatures are for SUSY models beyond 
the MSSM. In this work we consider SUSY left-right (LR) models with gauge mediated 
supersymmetry breaking (GMSB). 

Supersymmetric left-right models (SUSYLR) where the SU(2) R gauge symmetry is bro- 
ken by triplet Higgs fields A c with B — L = 2 have many attractive features: 1) they imply 
automatic conservation of baryon and lepton number 0; 2) they provide a natural solution 
to the strong and weak CP problems of the MSSM H; 3) they yield a natural embedding 
of the see-saw mechanism for small neutrino masses [|J where the right-handed triplet field 
(A c ) that breaks the SU(2)r symmetry also gives heavy mass to the right-handed Majorana 
neutrino needed for implementing the see-saw mechanism. 

Recently it has been shown that the doubly charged components of the triplet Higgs 
fields are massless unless there are some higher dimensional operators (HDO) ||5|-^J. This 
is independent of how supersymmetry breaking is transmitted to the visible sector (i.e., 
whether it is gravity mediated or it is gauge mediated) and also of whether the hidden 
sector supersymmetry breaking scale is above or below the Wr scale. In the presence of 
HDO's, they acquire masses of order ~ v R /M P \. Since the measurement of the Z-width at 
LEP and SLC implies that such particles must have a mass of at least 45 GeV, this puts 
a lower limit on the Wr scale of about 10 10 GeV or so. For Wr near this lower limit, the 
masses of the doubly charged particles are in the 100 GeV range. The rest of the particle 
spectrum below the Wr scale can be the same as that of the MSSM with a massive neutrino 
or it can have an extra pair of Higgs doublets in the 10 TeV range depending on the structure 
of the model. 

The ordering of the sparticles masses in the SUSYLR model makes the SUSY signature 
distinctive from other SUSY models. Of primary importance to the signal is the identity 
of the sparticle which is the next to lightest SUSY particle (NLSP). Here the NLSP can 
be the lightest neutralino, the lighter stau, or the light doubly charged Higgsino (which we 
henceforth call the deltino). Mass spectra for this type of model have been studied previously 
||. It was found that one of the Higgs fields has a coupling with the third generation 
charged leptons which reduces the third generation charged slepton masses. Because of this, 
in SUSYLR models with GMSB the lighter stau is predominantly the NLSP whenever the 
deltino is too massive to play that role. As a result, the decay chains of the SUSY particles 
typically lead to the lighter stau. The fx then decays into a r lepton and a gravitino (G) which 
escapes the detector undetected (leading to missing energy). Since the gravitino mass is on 
the order of eV, the emitted r will have high px enhancing its detection possibility. Moreover, 
pair production of the light doubly charged Higgsinos always produces four r leptons. When 
the fi is the NLSP, this occurs through A c±± — > f^r 1 * 1 followed by f\ — > r G. When the 
A c=t± is the NLSP, this occurs through the stau mediated decay A c±± — > r ± r ± G. One can 
get a similar signal in supergravity motivated LR models with the gravitinos replaced by 
the lightest neutralino (its mass is greater than 33 GeV), which will constitute the missing 
energy. 



2 



Signals involving two or more high t leptons are also important signals for conventional 
GMSB models as the lighter stau is frequently the NLSP for these models as well. In the 
SUSYLR model, however, we find that the production of the deltino can greatly enhance 
the signal. In addition, since the deltino decays into like sign r leptons, we find that the 
distribution in angle between same sign r leptons can be used to distinguish this model from 
other GMSB models. 

In this paper we analyze in detail the signal for SUSY production in the context of the 
SUSYLR model. First, we discuss the sparticle masses and the production cross sections 
for the major SUSY production modes: X1X1 production, xt X 2 production, slepton pair 
production and A c++ A c production. Next, we analyze the r-jet signal. We give the typical 
two, three and four r-jet inclusive production cross sections with and without standard cuts. 
Last, we show how the distribution in angle between the r-jets based on their charges can 
distinguish this model from other models whose signals also involve the production of r 
leptons. 



Since the observed signal depends on the masses of the sparticles, we first begin by 
describing the model and the corresponding mass spectrum. The particle content of this 
model above the LR scale includes 0(2, 2, 0), A(3, 1, 2), A(3, 1,-2), A c (l, 3, -2), A c (l, 3, -2) 
and a singlet where the numbers in the parentheses refer to their transformation properties 
under SU(2)l x SU(2)r x U(1)b-l- The LR symmetric superpotential for this theory is 



where Wnr denotes nonrenormalizable terms arising from higher scale physics such as grand 
unified theories or Planck scale effects. Demanding that the Fa, F s and F&c terms vanish, 
it has been found || that the VEV of S becomes ^y-. The VEVs for the A and A vanish 
and their masses become of the order of 2Mb- Let us now see what happens to the different 
components of the A c fields. One linear combination of the singly charged and neutral Higgs 
fields pick up a mass on the order of vr. The doubly charged Higgsinos are massless unless the 
nonrenormalizable terms are introduced. In the presence of Wnr, the doubly charged fields 
acquire masses of order v R /Mp\. This result is true when the vacuum conserves R-parity. 
When the vacuum breaks R-parity, the scale of the right handed PU's has an upper bound 
on the order of a few TeV. Therefore, in this version of the model, all particles, including 
the doubly charged bosons and fermions, have masses also in the few hundred GeV to TeV 
range. In what follows we work with the R-parity conserving version for simplicity. 

After integrating out the fields at the left-right scale, we are left with the following 
additional part to the MSSM: 



II. SPARTICLE MASSES AND PRODUCTION 



W = hfQ T T 2 ^r 2 Q c + h ( t ) L T T^ i T 2 L c 
+ i{iL T r 2 AL + f c L cT r 2 A c L c ) 

+ M D [Tr(AA) + Tr(A c A c )] + \S(AA - A C A T ) + ^ S S 2 

+ H ij TY(T2® T i T 2 ®j) + W NR 



(1) 



W = M A A C A c++ + fil c l c A c 



(2) 
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where we have assumed that f is diagonal. The PSI esperiment |T(J has put an upper bound 
on the product of the first two generation couplings of fif 2 < 1.2 x 1CT 3 . The magnitude of 
fs is unrestrained. The term M A originates from the nonrenormalizable terms. 

The model considered here involves gauge mediated supersymmetry breaking. In GMSB 
type models the SUSY breaking is communicated to the observable sector by the SM gauge 
interactions |ll| . We choose GMSB since the lighter stau will then have a mass that is 



almost always below that of the lightest neutralino. The lighter stau then decays to a r 
lepton and a gravitino. Since the gravitino is very light, the r lepton will typically be very 
energetic. 

In the GMSB model, the sparticle spectrum depends on the following parameters: M, 
A, n, tan (3, fa, M A (M) and the sign of fi. M is the messenger scale. The parameter n 
is dictated by the choice of the vector-like messenger sector. In this calculation we will 
assume that each flavor in the messenger sector consists of a vector like isosinglet pair of 
fields (Q + Q) and a vector like weak isodoublet pair L + L. The definition of tan/3 is 
that tan/3 = v 2 /vi where v 2 is the VEV for the up-type (H u ) Higgs doublet and V\ is the 
VEV for the down-type (Hd) Higgs doublet. M A (M) is the messenger scale value for the 
deltino mass. The parameter /i is the coefficient in the bilinear mixing term, /j,H u Hd, in the 
superpotential. Constraints coming from b — ► 57 strongly favor negative values for fi [T2 



and, in the cases considered in this work, /1 is taken to be negative. Demanding that the 
EW symmetry be broken radiatively fixes the magnitude of fi and the parameter B (from 
the BnH u Hd term in the scalar potential) in terms of the other parameters of the theory. 
The soft SUSY breaking gaugino and scalar masses at the messenger scale are given by 

H 

M i( M)= n9 fA^ A (3) 

and 

m 2 



MJ 4vr 



(M)^/(A)g,^)V (4) 



where the a>i are the three SM gauge couplings and fcj = 1, 1 and 3/5 for SU(3), SU(2) 
and U(l), respectively. The Ci are zero for gauge singlets and are 4/3, 3/4 and (Y/2) 2 for 
the fundamental representations of SU(3), SU(2) and U(l), respectively (with Y given by 
Q = I 3 + Y/2). g(x) and f(x) are messenger scale threshold functions. 

We calculate the SUSY mass spectrum by using the appropriate renormalization group 
equations jl4]| . We first run the Yukawa couplings (including the three new couplings /i,2,3) 



and the gauge couplings from the weak scale up to the messenger scale. At the messenger 
scale, we apply the boundary conditions given by the equations above and then use the 
RGEs for the soft SUSY couplings and masses in order to run down to the weak scale. 

The mass spectrum here is much like that expected in minimal GMSB models. The 
gravitino is always the LSP. Since SUSY breaking is communicated to the visible sector by 
gauge interactions, the mass differences between the superparticles depend on their gauge 
interactions. This creates a hierarchy in mass between electroweak and strongly interacting 
sparticles. Eq. |3| shows that the gluino is more massive than the charginos and neutralinos, 
while Eq. |] shows that the squarks are considerably more massive than the sleptons. Thus 
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in minimal GMSB models, the lightest neutralino and the lighter stau fight for the NLSP 



spot [|12yi5l . In this model, the deltino also joins the race to become the NLSP. 

We will concentrate the analysis on those regions of the parameter space where either 
the lighter stau or the deltino is the NLSP. Whether or not the deltino is the NLSP depends 
on the mass it gets from the higher dimensional terms. If this mass is too high, then either 
the fi or Xi is the NLSP. The lighter stau can be much lighter in our SUSYLR model than 
in conventional GMSB models due to the presence of the additional coupling ^3. Thus the 
lighter stau will be lighter than the Xx f° r a larger region of the parameter space and the fx 
has a greater potential to be the NLSP in this SUSYLR model. 

There are a number of potential SUSY production mechanisms here. Given the current 
lower bounds on the various sparticle masses and the hierarchy of sparticle masses in GMSB 
models, the important SUSY production mechanisms will typically include EW gaugino 
production. At the Tevatron, chargino pair (xxXi) production takes place through s-channel 
Z and 7 exchange and X2 Xx production is through s-channel W exchange. Squark exchange 
via the t-channel also contributes to both processes, but the contributions are expected to 
be negligible since the squark masses are large in GMSB models. The production of Xi Xi 
is suppressed due to the smallness of the coupling involved. 

In addition to these usual SUSY production mechanisms of the MSSM, we also have 
deltino pair (A c++ A c ) production. This proceeds through s-channel Z and 7 exchange. 
Given that the A c±± can be relatively light, it can be a very important SUSY production 
mode. In fact, it frequently is the dominant mode. 

The possible final state configurations at the Tevatron depend on the sparticle spectrum 
and on which SUSY production mode is dominant, but they will have certain aspects in 
common. When the fx is the NLSP, the various possible decay modes will (usually) produce 
at least two r leptons arising from the decays of the lighter staus. In addition, there can also 
be large JpT due to the stable gravitinos and neutrinos escaping detection. When the deltino 
is the NLSP, the standard SUSY production modes of EW gauginos can still produce large 
numbers of r-jets if the f\ is the next to next to lightest SUSY particle (in which case the 
fx is lighter than the Xx) so that the decay chains of the sparticles will still lead to the f\. 

Pair production of the deltino leads to copious quantities of r leptons irrespective of 
what the NLSP is. This is because the deltino couples to leptons/sleptons but not to 
quarks/ squarks. In addition, the coupling to the third generation can be much greater 
than the small coupling to the 1st and 2nd generations. Thus, when the fx is the NLSP, 
the deltino decays via A c±± — > f^r^ with the stau decaying via fi — > tG. On the other 
hand, when the deltino is the NLSP, it decays via the f mediated three-body decay mode 
A c±± — > r ± r ± G'. Thus A c±± pair production generally results in the production of four r 
leptons (two from each deltino). 



III. TAU JET ANALYSIS 



As mentioned above, SUSY production for this SUSYLR model leads to the production 
of copious quantities of r leptons. r leptons are typically identified at colliders by their 
hadronic decays to thin jets. We now give a detailed account of the possible r-jet signatures 
for SUSY production at the Tevatron in the context of the left-right GMSB model. 
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This analysis is performed in the context of the Main Injector (MI) and TeV33 upgrades 
of the Tevatron collider. The center of mass energy is taken to be y^i = 2TeV and the 
integrated luminosity is taken to be 2fb _1 for the MI upgrade and 30 fb -1 for the TeV33 
upgrade. 

In performing this analysis, the cuts employed are that final state charged leptons must 
have pt > 10 GeV and a pseudorapidity, r\ = — ln(tan |) where 6 is the polar angle with 
respect to the proton beam direction, of magnitude less than 1. Jets must have Et > 10 GeV 
and < 2. In addition, hadronic final states within a cone size of AR = \J (A0) 2 + (Arj) 2 = 
0.4 are merged to a single jet. Leptons within this cone radius of a jet are discounted. For 
a r-jet to be counted as such, it must have \r]\ < 1. The most energetic r-jet is required to 
have Et > 20 GeV. In addition, a missing transverse energy cut of $t > 30 GeV is imposed. 

The signatures for SUSY production depend on the hierarchy of sparticle masses. This, in 
turn, depends on the values the parameters of the theory takes. The parameters considered 
in this analysis are tan (3 = 15, n = 2, M/A = 3, / 3 = 0.5, f 2 = 0.05 and f x = 0.05. We 
vary A from 35 to 85 TeV. For the messenger scale deltino mass, we use the values 90, 120 
and 150 GeV. The masses of some of the particles of interest are given in Fig. I and Fig. |. 
In Fig. [I] we take M^(M) = 90 GeV, but the masses of the gauginos and sleptons (with the 
exception of the stau) do not vary much with the messenger scale deltino mass. Fig. ||] gives 
the masses of the delta boson and the deltino. The deltino mass is not very sensitive to the 
value of A, while the delta boson mass is highly dependent on A due to the contributions 
from the messenger scale loops (which give it mass along with the nonrenormalizable part). 
Given the substantially higher A c boson mass, A c production is not very important at the 
Tevatron. 

There are several potential SUSY production modes here. The cross sections for the 
more traditional SUSY production modes are given in Fig. [| We also have deltino pair 
production; the cross sections for which are tabulated in Table [I| Since the deltino mass does 
not vary much over the values of A considered, the cross section for deltino pair production 
does not vary much either. This cross section is high enough for all the deltino masses 
considered that deltino pair production is always an important SUSY production mode. 
For low values of A, the EW gaugino production cross section is large with values in the 
hundreds of fb at A = 35 TeV, but the cross section falls off substantially as A increases. As 
A increases above about 55 TeV, the cross section for EW gaugino production starts to fall 
below that of slepton production (in particular f^rf ). In a minimal model, these sleptons 
modes would become the dominant SUSY production modes, but here the cross sections for 
slepton production fall far below that of deltino pair production. Thus the dominant SUSY 
production modes here are deltino pair production and, at values of A below 45 TeV or so, 
XiXi and \\ Xi production. 

The decay chains depend on which sparticle is the NLSP. For the values of the parameters 
that are considered here, either the lighter stau or the deltino is the NLSP. Since the mass 
of the lighter stau increases with increasing A, the lighter stau is the NLSP for lower values 
of A, while the deltino is the NLSP for higher values of A. For a messenger scale deltino 
mass of 90 GeV, the lighter stau is the NLSP for A below about 43 TeV. For M A (M) = 120 
and 150 GeV, the boundaries are given by about 58 and 73 TeV, respectively. When the 
ligher stau is the NLSP, it decays via fi — > tG, and the deltino decays via the two-body 
mode A c — > f\T. Then deltino pair production leads to four r leptons. On the other hand, 
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if the deltino is the NLSP, it decays via the stau mediated three-body mode A c — > ttG. So, 
once again, deltino pair production again leads to the production of four r leptons. 

The decays of the lighter selectron and smuon are given in Table O. At values of A 
around 35 to 40TeV, the neutralino is lower in mass than the e\ and fix. When this is 
the case, the main decay mode of the is e.\ — > Xi e an d the smuon decay is correspondingly 
fii Xi A*. When this decay is not kinematically allowed, the decay e.\ — > A c e is typically 
dominant if kinematically allowed. If it isn't, then the selectron decays via the three-body 
decays ef — > e + r + ff and ef — > e + T~?x and/or the two-body mode e.\ — > eG. 

The branching ratios for the neutralinos and lighter chargino are given in Table III. The 
lighter neutralino has only the three decay modes Xi ~^ ^i r > Xi A^i/ i an d Xi ~ > &i e over 
the parameter space considered. Since the lighter neutralino is lighter than fix and e.\ for 
A < 43 TeV, the only decay mode for Xi is Xi ~~ * ^i r - As A increases beyond the point 
where the decays to the selectron and smuon become kinematically available, the branching 
ratios for Xx ~^ fiiH an d X? — > &i e increase, but the Xi decay remains dominant due 

in large part to the fact that the mass of f\ is much lower than that of the selectron and 
smuon. 

The chargino has only two decay modes over the allowed parameter space: Xi ^i u t 
and Xi ~^ Xi W. At the lower values of A considered, the decay to the lighter stau is either 
the only decay mode available or is the dominant decay mode. For A around 40 TeV and 
below, the only decay mode for the lighter chargino is Xi ~^ ^\ v t- F° r these values of A, 
the lighter stau decays via f\ — > r G as discussed above. Thus in chargino pair production, 
two r leptons are produced. As A increases, the decay mode Xi ~ > Xi W appears. With 
the subsequent decays of the lighter neutralino to the sleptons and with the deltino as the 
NLSP, the number of r leptons produced is typically four or six (in principle eight r leptons 
can be produced although this requires the rather rare three-body decays of the selectron 
and smuon). 

The branching ratios of the second lightest neutralino are particularly sensitive to the 
value of A. At lower values of A, the decays to the sleptons are dominant. In particular, 
the decay X2 ~^ i s dominant due to the lower mass of the f\ and the fact that the x\ 
is mostly wino. When the t\ is the NLSP, x% Xi production typically produces three r-jets. 
When the lighter stau isn't the NLSP and decays via f\ — > r A c , then five r leptons are 
usually produced, although three is also common due to the decays of the X2 t° the fix and 
ei followed by their decays to the deltino. As A increases, the decay X2 ~* Xi h becomes 
dominant, but at these values of A, the cross section for EW gaugino production falls far 
below that of deltino pair production. 

In summary, the dominant SUSY production modes at low values of A are deltino pair 
production and EW gaugino production. We expect four r leptons to be produced in deltino 
pair production, while EW gaugino production is typically expected to produce two to three 
r leptons. For larger values of A, the possibility exists to produce many r leptons in EW 
gaugino production, but the cross sections for such production are much smaller than that 
for deltino pair production. Thus four r leptons are generally produced at larger values of 
A. 

We now consider the observability of these modes at Tevatron's Run II. Tables [V], [V] 
and [VI] give the inclusive r-jet production cross sections for a messenger scale deltino mass 
of 90, 120 and 150 GeV, respectively. We include in the figures only up to four r-jets as the 



7 



cross sections for more than four r-jets are small. Considering Table [IV], we see that before 
cuts the production of two and three r-jets are dominant, but the four r-jet cross section 
is also significant at slightly over 100 fb. After the cuts are applied, however, the situation 
changes substantially. The one r-jet mode is now dominant, but the cross section for two 
r-jets is not far below and the three r-jets cross section is not insignificant. 

We first consider the M^(M) = 90GeV case. We see from Table [IV] that for A = 
35TeV the cross section for inclusive production of three r-jets is 32.3 fb. For an integrated 
luminosity of 2fb _1 (the approximate intial value at Run II), this corresponds to about 
65 events. For 30 fb -1 , the number of observable events is ~ 970. For A = 85TeV, the 
production cross section for three r-jets has gone down slightly due to the decrease in 
production of charginos and neutralinos. With a value of 26.7 fb, the number of expected 
events is about 53 and 800 for 2fb _1 and 30 fb -1 of data, respectively. The cross section 
for two r-jets is considerably higher. For A = 35 TeV, the a • BR for two r-jets is 125 fb 
which corresponds to 250 events for 2fb _1 of data and 3750 events for 30 fb -1 of data. 
For A = 85 TeV, a ■ BR has decreased to 79 fb. This gives about 160 and 2370 events for 
2fb _1 and 30 fb -1 of data, respectively. In comparison to the GMSB model with the MSSM 
symmetry, the two r-jets and the three r-jets cross sections are considerably higher in this 
model. In the GMSB model with MSSM symmetry, the two r-jets cross section can be seen 



at RUN II, but not the three r jets [16 



As the mass of the deltino increases, the production rates go down and more variations 
appear. The inclusive r-jet cross sections for M^(M) = 120 GeV are shown in Table [V]. 
Considering the inclusive three r-jets mode, the production cross section at A = 35 TeV is 
26 fb. This corresponds to 52 events for 2ib~ 1 of data and 780 events for 30 fb -1 of data. 
This goes down to about 17 fb at A = 85 TeV. This gives about 34 and 510 events for 
2fb _1 and 30 fb -1 of data, respectively. The production rate for two r-jets is higher. At 
A = 35 TeV, a ■ BR = 94 fb which gives 190 and 2820 events for 21b" 1 and 30 fb" 1 of data, 
respectively. 



IV. ANGULAR DISTRIBUTIONS 

The excess of r-jets expected in this model does not constitute an unequivocal signal for 
this model, r-jets are part of the signatures for other models including the minimal GMSB 
model when the lighter stau is the NLSP. The question then arises as to whether there is 
any way to distinguish this model from the minimal GMSB model. A possible distinguishing 
characteristic is the distribution in angle between the two highest r-jets when they come 
from same sign r-jets. 

Consider deltino pair production. The deltino tends to decay to like sign r leptons. 
This occurs directly when the deltino is the NLSP and so decays via the three-body decay 
— > r ± r ± G'. When the two-body decay of the deltino A ±=l= — > r^r 1 * 1 occurs, then the 
second like sign r lepton comes from the subsequent decay of the stau. In the rest frame 
of the deltino, the r leptons are widely distributed. In the lab frame, however, the deltinos 
are quite energetic and have a large velocity, especially if their masses are small. As a 
consequence of this, the decay products of the deltino tend to be collimated in the direction 
in which the deltino was moving. Thus when the two most energetic r-jets have the same 
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sign in deltino pair production, the angle between them tends to be smaller than when the 
two most energetic r-jets have opposite sign charges. 

Fig. gives the distribution in angle between the two most energetic r-jets for deltino 
pair production. This example is for a weak scale deltino mass of about 97GeV. We can 
see that the distribution in angle for like sign r-jets, which is given in Fig. f|(a), peaks at 
about 40°. Fig. f|(b) gives the distribution in angle between the two most energetic r-jets 
when they come from opposite sign r leptons. In stark contrast to the previous case, here 
the peak occurs at 110°. 

The question then arises as to how these distributions look in the usual SUSY production 
modes. Fig. |5] shows the angular distributions for combined X2 Xi an d xtxl production for 
the input parameters M = 100 TeV, A = 45 TeV, n = 1 and tan (3 = 10. For these values 
of the parameters, the weak scale X2 m &ss is ~ 100 GeV. The distribution for same sign 
r-jets is given in Fig ||(a). We see that the peak occurs at about 110°. In this situation, 
same-sign r-jets do not come from X1X1 production. In \\ Xi production, one of the same 
sign r-jets generally comes from the chargino and the other from the neutralino. We now 
consider the angular distribution for opposite sign r-jets which are given in Fig. |^(b). In 
X2 Xi production, opposite sign r-jets frequently come from the neutralino, while in X1X1 
production one of the r-jets comes from one of the charginos and the other r-jet comes from 
the other chargino. Since there is a strong possibility that the opposite sign r-jets come 
from the same particle {x^i the distribution should peak at a lower angle than for same 
sign r-jets. We see from the figure that the peak occurs at about 85°. 

The question arises as to how much this changes as the gaugino masses are increased. 
Fig. |6| gives the angular distribution for a x\ mass of 150 GeV. We see that the same sign 
distribution still peaks at about 110°, while the opposite sign distribution has now shifted 
to a slightly higher value of about 95°. 

The actual angular distribution between the two highest Et r-jets depends on which 
SUSY production modes are important. For certain regions of the parameter space (de- 
pending, in particular, on the values of A and the messenger scale deltino mass), deltino 
pair production is the only important SUSY production mode. When this is the case, the 
angular distributions are simply given by those for deltino pair production. 

In other regions of the parameter space, EW gaugino production can be significant. We 
consider the angular distributions for an example given by the input parameters tan/3 = 15, 
n = 2 and M/A = 3. Three values of the messenger scale deltino mass are considered: 90, 
120 and 150 GeV. The angular distributions for M-^(M) = 90 GeV are given in Fig. [7|. Since 
the deltino is especially light at 96 GeV, deltino pair production is the dominant SUSY 
production mode. Thus deltino pair production largely dictates the form of the angular 
distributions. 

Fig. 0(a) gives the angular distribution between the two highest Et r-jets when they 
come from same sign r leptons. In the figure we can see the rather striking peak at around 
40°. This is due to the same sign r-jets coming mostly from the decay of the same deltino. 
Since the deltino mass is especially light compared to the beam energy, they typically move 
rapidly in the lab frame. Thus their decay products tend to be more tightly collimated than 
in the production of the heavier particles. 

Fig. 0(b) gives the angular distribution between the two highest Et r-jets when they 
come from opposite sign r leptons. We see from this figure that the peak occurs at about 
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110°. Here the r-jets typically come from the decay chains of different particles and so the 
angle between the r-jets is typically quite large. 

The situation changes as the deltino mass gets larger. This is in part due to the fact 
that the deltino pair production cross section gets smaller and so production of charginos 
and neutralinos can have a larger impact on the distributions. In addition, a larger deltino 
mass means the deltinos will typically be moving slower. Thus the boost effect won't have as 
dramatic an effect on the deltino's decay products. The example with M^(M) = 120 GeV 
is given in Fig. §. We can see that the distribution for same sign r-jets peaks at about 70°. 
On the other hand, the opposite sign r-jet angular distribution still peaks at around 110°. 
Thus the angle between the r-jets is less striking a signature than it was before, but it is 
still distinctive. 

The results for a messenger scale deltino mass of 150 GeV is given in Fig. |9|. The peak in 
the distribution in angle between the two highest r-jets when they have the same sign 
peaks at a rather high 95°. As before, the peak in the distribution for the two highest E? 
r-jets when they have opposite sign is at 110°. Thus the distinctiveness due to the angle 
between the two highest Et r-jets is nearly lost for such a large value of the deltino mass. 
This is partially due to the deltino pair cross section of 92 fb being quite a bit lower than the 
cross section for Xi X2 production and X1X1 production. In addition, as discussed above, 
there is also the reduction in the boost effect as the mass of the decaying deltino increases. 

V. CONCLUSION 

In conclusion, we have found that the doubly charged Higgs bosons of LR models can 
be potentially observable at Run II of the Tevatron through the production of r-jets. In 
a GMSB type theory, SUSYLR models typically produce large numbers of two and three 
r-jet final states. This large r-jet signal is also due in large part to pair production of the 
doubly charged Higgsino. It is also due to the relatively low mass of the lighter stau (which 
is frequently the NLSP) in these models, which is due to the additional coupling /. We have 
also shown that the distribution in angle between the two highest E? r-jets is different from 
other models which do not have this doubly charged Higgsino. 
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FIGURES 
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FIG. 1. Masses of the particles of interest for the input parameters tan/3 = 15, M/A = 3, 
n = 2, / 3 = 0.5, h = 0-05, fi = 0.05 and M A (M) = 90GeV. 
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FIG. 2. Masses of the delta and deltino. The dashed lines represent the deltino, while the 
solid lines represent the delta boson. The parameters used are tan/3 = 15, n = 2 and M/A = 3. 
From bottom to top, the lines in each set are for a messenger scale deltino mass of 90, 120 and 
150 GeV. 
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FIG. 3. Cross sections for the standard SUSY production modes for the parameters tan /3 = 15, 
M/A = 3, n = 2, f 3 = 0.5, f 2 = 0.05, /i = 0.05 and M K {M) = 90 GeV. 
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TABLES 



TABLE I. Cross sections (in fb) for deltino pair production for various values of A. The other 
parameters used are tan (3 = 15, n = 2 and M/A = 3. 



M A (M) 




A = 35 TeV 


A: 


= 60 TeV 




A = 85 TeV 


90 GeV 
120 GeV 
150 GeV 




643.0 
228.2 
91.6 




629.9 
222.9 
89.2 




621.5 
219.6 
87.7 


TABLE II. 

n = 2 and M/A 


Branching ratios of the sleptons. The values of the parameters 
= 3. 


are tan/3 = 15, 


A (TeV) 


35 


40 


50 


60 


70 


80 85 


M A (M) = 90 GeV 


ei -> ex? 
ei — > e A 


1 


0.9252 
0.0748 


1 


1 


1 


1 1 


fl -4TG 

fi — ► r A 


1 


1 


1 


1 


1 


1 1 


M A (M) = 120 GeV 


ei -> ex? 
— > e + r + ff 
— > e+r^f^ 

ei — > e A 


1 


1 


0.2065 
0.1684 
0.6251 


1 


1 


1 1 


n^rG 
fi — > r A 


1 


1 


1 


1 


1 


1 1 


M A (M) = 150 GeV 


ei -> ex? 
— > e + T + ff 
— > e + T~fi~ 

ei — > e A 


1 


1 


0.5532 
0.4468 


0.5643 
0.4357 


1 


1 1 


h^rG 
fi — > r A 


1 


1 


1 


1 


1 


1 1 



15 



TABLE III. Branching ratios of some of the sparticles of interest. The values of the parameters 
are tan (3 = 15, n = 2 and M/A = 3. The messenger scale deltino mass is 90 GeV, but the branching 
ratios of these sparticles have little dependence on the deltino mass. 



A (TeV) 


35 


40 


50 


60 


70 


80 


85 




1 


1 


0.7153 


0.5807 


0.5166 


0.4796 


0.4663 


xt Xi w 






0.2847 


0.4193 


0.4834 


0.5204 


0.5337 


X° hr 


0.6312 


0.6587 


0.6860 


0.4034 


0.3133 


0.2706 


0.2560 


X° g ie 


0.1844 


0.1707 


0.1366 


0.0625 


0.0386 


0.0270 


0.0232 


X° AiM 


0.1844 


0.1707 


0.1366 


0.0625 


0.0386 


0.0270 


0.0232 


xS x? z 






0.0408 


0.0325 


0.0277 


0.0252 


0.0243 


X2 Xi h 








0.4392 


0.5818 


0.6501 


0.6733 


Xi TiT 


1 


1 


0.9692 


0.9099 


0.8723 


0.8506 


0.8436 


Xi eie 






0.0154 


0.0451 


0.0638 


0.0747 


0.0782 


Xi AiM 






0.0154 


0.0451 


0.0638 


0.0747 


0.0782 



TABLE IV. Inclusive r-jet production cross sections for a messenger scale deltino mass of 
90 GeV. The other parameters are tan/3 = 15, n = 2 and M/A = 3. 



A (TeV) 


35 


40 


50 


60 


70 


80 


85 










a ■ BR (fb) 








1 r-jet: before cuts 


198.6 


132.9 


73.11 


71.64 


70.86 


70.34 


70.14 


after cuts 


168.7 


156.2 


97.49 


91.90 


90.15 


89.33 


89.05 


2 r-jets: before cuts 


362.6 


277.3 


203.6 


196.7 


196.3 


194.8 


194.2 


after cuts 


124.6 


93.11 


89.59 


82.91 


80.59 


79.63 


79.33 


3 r-jets: before cuts 


306.0 


273.5 


253.6 


244.7 


240.9 


238.8 


238.0 


after cuts 


32.31 


17.75 


32.45 


29.02 


27.53 


26.88 


26.73 


4 r-jets: before cuts 


119.4 


116.9 


126.42 


116.5 


113.1 


111.6 


111.1 


after cuts 


3.21 


1.18 


5.26 


4.28 


3.67 


3.45 


3.39 
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TABLE V. Inclusive r-jet production cross sections for a messenger scale deltino mass of 
120 GeV. The other parameters are tan/3 = 15, n = 2 and M/A = 3. 



A (lev ) 




/in 


ou 


ou 


I U 


oU 


QK 










a ■ Dtt (ID J 








1 r-jet: before cuts 


151.5 


86.21 


39.60 


29.72 


25.30 


25.00 


24.90 


after cuts 


125.3 


85.20 


53.68 


75.87 


42.16 


41.38 


41.16 


2 r-jets: before cuts 


232.8 


148.0 


90.27 


75.93 


70.06 


69.15 


68.83 


after cuts 


93.94 


71.59 


49.79 


28.57 


43.81 


42.83 


42.55 


3 r-jets: before cuts 


147.6 


115.9 


96.70 


88.92 


86.55 


85.16 


84.71 


after cuts 


25.99 


23.22 


17.33 


0.55 


17.37 


16.80 


16.64 


4 r-jets: before cuts 


46.22 


44.03 


42.59 


40.70 


40.72 


39.55 


39.22 


after cuts 


3.05 


3.09 


2.29 


0.06 


2.70 


2.50 


2.44 



TABLE VI. Inclusive r-jet production cross sections for a messenger scale deltino mass of 
150 GeV. The other parameters are tan/3 = 15, n = 2 and M/A = 3. 



A (TeV) 


35 


40 


50 


60 


70 


80 


85 










a • BR (fb) 








1 r-jet: before cuts 


136.2 


70.92 


24.42 


14.61 


11.80 


10.07 


10.00 


after cuts 


105.3 


65.54 


32.12 


22.59 


23.50 


19.28 


18.97 


2 r-jets: before cuts 


190.3 


105.7 


48.32 


34.22 


30.01 


27.89 


27.67 


after cuts 


72.17 


50.55 


31.35 


24.29 


18.99 


21.75 


21.46 


3 r-jets: before cuts 


95.31 


63.85 


45.21 


37.72 


35.14 


34.47 


34.14 


after cuts 


17.23 


14.64 


12.24 


10.14 


5.53 


9.32 


9.24 


4 r-jets: before cuts 


21.81 


19.75 


18.58 


17.00 


16.16 


16.23 


15.96 


after cuts 


1.73 


1.81 


1.89 


1.67 


0.59 


1.46 


1.47 
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0.012 




FIG. 4. Angular distribution between the two most energetic r-jets for deltino pair production 
at the Tevatron. The deltino mass is about 97 GeV. (a) gives the distribution when the r-jets come 
from same sign r leptons. (b) gives the distribution when the r-jets come from opposite sign r 
leptons. 
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FIG. 5. Angular distribution between the two most energetic r-jets for EW gaugino production 
at the Tevatron. The X2 mass is 100 GeV. (a) gives the distribution when the r-jets come from 
same sign r leptons. (b) gives the distribution when the r-jets come from opposite sign r leptons. 
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FIG. 6. Angular distribution between the two most energetic r-jets for EW gaugino production 
at the Tevatron. The mass of \2 IS about 150 GeV. (a) gives the distribution when the r-jets come 
from same sign r leptons. (b) gives the distribution when the r-jets come from opposite sign r 
leptons. 
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FIG. 7. Angular distribution between the two most energetic r-jets for combined SUSY pair 
production at the Tevatron. The messenger scale deltino mass is 90 GeV. The other parameters 
are tan/3 = 15, n = 2 and M/A = 3. (a) gives the distribution when the r-jets come from same 
sign r leptons. (b) gives the distribution when the r-jets come from opposite sign r leptons. 
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FIG. 8. Angular distribution between the two most energetic r-jets for combined SUSY 
production at the Tevatron. The messenger scale deltino mass is 120 GeV. The other parameters 
are tan/3 = 15, n = 2 and M/A = 3. (a) gives the distribution when the r-jets come from same 
sign r leptons. (b) gives the distribution when the r-jets come from opposite sign r leptons. 
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FIG. 9. Angular distribution between the two most energetic r-jets for combined SUSY 
production at the Tevatron. The messenger scale deltino mass is 150 GeV. The other parameters 
are tan/3 = 15, n = 2 and M/A = 3. (a) gives the distribution when the r-jets come from same 
sign r leptons. (b) gives the distribution when the r-jets come from opposite sign r leptons. 
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